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Abstract 

In this paper we study the integrability of a class of Gross-Pitaevskii equations managed 
by Feshbach resonance in an expulsive parabolic external potential. By using WTC test, 
we find a condition under which the Gross-Pitaevskii equation is completely integrable. 
Under the present model, this integrability condition is completely consistent with that pro- 
posed by Serkin, Hasegawa, and Belyaeva [V. N. Serkin et al., Phys. Rev. Lett. 98, 074102 
(2007)]. Furthermore, this integrability can also be explicitly shown by a transformation, 
which can convert the Gross-Pitaevskii equation into the well-known standard nonlinear 
Schrodinger equation. By this transformation, each exact solution of the standard nonlinear 
Schrodinger equation can be converted into that of the Gross-Pitaevskii equation, which 
builds a systematical connection between the canonical solitons and the so-called nonau- 
tonomous ones. The finding of this transformation has a significant contribution to under- 
standing the essential properties of the nonautonomous solitions and the dynamics of the 
Bose-Einstein condensates by using the Feshbach resonance technique. 
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1 Introduction 



Dilute-gas Bose-Einstein condensates (BECs) have been generated by many exper- 
iment groups since 1995 (1; 2). It is well-known that the mean-field dynamics of 
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the BECs at low temperature can be well described by a three-dimensional Gross- 
Pitaevskii (GP) equation (3; 4), where the BECs are confined to an expulsive har- 
monic external potential. In some physically important cases the GP equation can 
be reduced effectively to a one-dimensional GP equation (5; 6; 7; 8; 9; 10; 11; 12). 
This simplified equation reads after re-scaling the physical variables, 



where i(j(x, t) denotes the macroscopic wave function of the condensate and V(x) is 
the potential confining the condensate. The g represents the nonlinear interaction 
strength and a time-dependent factor of T(t) is controlled by the so-called Fesh- 
bach resonance(13; 14). It is very interesting to study the dynamics of the BECs 
described by Eq. (1) with the Feshbach resonance management. 

It should be pointed out that Eq. (1) is a special case of a more generalized nonlinear 
Schrodinger (NLS) equation with varying in time and space dispersion and nonlin- 
earity. In the literature, this generalized NLS equation was named as the nonau- 
tonomous NLS equation (15). The exact solutions of the nonautonomous NLS 
equation are suggested to name the nonautonomous solitons in order to conception- 
ally distinguish from the canonical solitions introduced by Zabusky and Kruskal 
(16). Actually, the studies of the integrability of the generalized NLS equation have 
a long history. As early as in 1976 Chen and Liu (17) went beyond the concept 
of the canonical soliton and found the exact integrability and the Lax pair of NLS 
equation with space-varying potentials. More recently, the dynamics of the solitons 
of these generalized NLS equation has been extensively studied (15; 18; 19) and 
some novel concepts, i.e., the nonautomonous NLS equation (15) and the nonau- 
tomonous soliton (19) have been proposed. 

A classical way to study the integrability of the standard NLS equation and the 
nonautomonous one is the construction of the Lax pair of these NLS equations 
based on the inverse scattering transform (20). In this paper, we use the WTC test 
method suggested by Weiss, Tabor and Carnevale (21; 22) in 1983 to study the in- 
tegrability of the GP equation with Feshbach resonance management. This method 
is based on the Painleve test for PDE. Through the WTC test, we get an integrabil- 
ity condition. At the same time, we also find a transformation that can convert the 
GP equation with Feshbach resonance management into the standard NLS equa- 
tion, which has been extensively studied and a number of exact solutions including 
various solitons have been obtained in the literature. From the well-known solu- 
tions of the standard NLS equation, we can obtain novel solutions of the integrable 
GP equation with Feshbach resonance management. This transformation builds a 
systematical connection between the well-known canonical soltions and the nonau- 
tomonous ones and provide a novel way to study the dynamics of the BEC de- 
scribed by the GP equation with Feshbach resonance management. Finally, we also 
find a new type of solution of the standard NLS equation, which is also a new type 




(1) 



2 



of solution of the GP equation. 



The outline of this paper is as follows. In Sec. 2 we explicitly present the WTC test 
of the GP equation (1) and discuss the condition the equation can pass the WTC 
test, i.e., the integrability condition of Eq. (1). In Sec. 3 we show a transformation 
to convert the GP equation (1) into the well-known standard NLS equation, which 
demonstrates the complete integrability of the GP equation under the condition 
obtained in the previous section. In Sec. 4, as some applications, we list some solu- 
tions of the standard NLS equation, including a new type of solution, from which 
numerous novel solutions of the integrable GP equation with Feshbach resonance 
management can be obtained by the transformation. 



2 The Painleve test for the GP equation with Feshbach resonance manage- 
ment 



For simplicity, we rewrite Eq. (1) as 

d d 2 

i—u (x, t) + — — -u (x, t) + g (x, t) \u (x, t) | u (x, t) + V x 2 u (x, t) — (2) 



dt v 1 ' dx 2 

In order to use the Painleve analysis, it is convenient to introduce a complex func- 
tion v(x, t) = u(x, t)* . Thus, Eq.(2) becomes a pair of complex equations, 

ij^u (x, t) + -j^u (x, t) + g (x, t) (u (x, t)) 2 v (x, t) + Vqx 2 u (x, t) = 0, 
-i-^v {x, t) + (x, t) + h 0, t) (v (x, t)) 2 u (x, t) + V x 2 v (x, t) = 0, 

where h(x, t) = g(x, t)* . 

The next step is to seek a solution of Eqs.(3) in the following form 

u (x, *) = 4>(x, t)~p EZo ^ (*) (4 

v (x, t) = (f>( x , tyi et=o Vj{t)<t>{x, ty. 

The coefficients g(x, t) and h(x, t) can be expanded on the singularity manifold as 
follows 

' g(x,t) = Y,%Lo9k{t)<f>{x,t) k , 

h(x,t) = Zn l=0 h m {t)<P{x,t) m . 

Here we use the Kruskal assumption (f)(x, t) = x+ip(t), where Ui(t), Vj(t), gk(t), h m (t) 
are analytic functions in the neighborhood of a non-characteristic singularity man- 
ifold defined by x = — ip(t). 
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By the standard procedure, the leading order is given by p = q = 
obtain the following relations 



1, and one can 



Qi 9o (t) (u (t) f 




Uj (t) 






h (t)(v (t)) 2 Q 2 








G 3 _ 



where j = 1, 2, • • • , and 



(6) 



Qi = {j - 1) (j - 2) + 2 g (t) u (t) v (t) , 
Q2 = {j - 1) {j -2) +2 h (t) uo (t) u (t) • 

Here #o(£)> ^o(*)> «o(*) and v (t) obey 

#0 (*) % (t) v (t) = h (t) u (t) v (t) = -2. 
The recursion relations are determined by 



(7) 



Fj = -jUj-2 (t) - iuj-i (t) (j^ (t)j j + 2 iuj-i (t) (t) - 9j (t) (u (t)) 2 v (t) 

j—l j-l m 

~90 (t) V (t) Y U j-m (t) U m (t) - g (t) Y V j-m if) Y Um-k (t) U k (t) 
m=l m=l k=0 

j—l m I 

~ 9j-m (t) Y V m-l (t) Y U l~k (*) U k (*) ~ V X 2 Uj- 2 (*) j (8) 
m=l Z=0 fc=0 

Gj = ~j t Vj-2 (t) + iVj-! (t) (t)\ j - 2 iv^ (t) (t) - hj (t) (v (t)) 2 Uo it) 

j—l j-l m 

-h (t) U (t) Y V 3-m (t) V m (t) - h (t) Y Uj- m (t) Y V m -k (t) V k (t) 
m=l m=l k=0 

j—l m I 

- hj-m (t) Y U m-l (t) Y V l-k (t) ^ (t) ~ V X 2 Vj - 2 (t) • (9) 

m=l 1=0 k=0 

One notes that gj(t) = if j < 0. It is also true for hj(t),Uj(t), Vj(t). The expres- 
sions Fj, Gj for a given j depend only on the expansion coefficients of u\, vi, with 
I < j and gi, hi with I < j. Therefore the above equations represent the recursion 
relations of the unknown Uj, Vj (j > I) from the known ui, vi (I < j). 

The above recursion relations determine the unknown expansion coefficients uniquely 
unless the determinant of the matrix in Eq. (6) is zero. Those values of j at which 
the determinant is equal to zero are called as the resonance points. After some cal- 
culations it is found that the resonance points only occur at 



J 



-1,0,3,4. 
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The resonance point of j = — 1 corresponds to the arbitrariness of the singular 
manifold (f)(x, t). 

From the recursion relations, we find the compatibility conditions for the remaining 
resonance points 



3=0: g (t) u (t) v (t) = h (t) u (t) v (t) = -2, (10) 
j = 3: v (t)F 3 = u (t)G 3 , (11) 
3 = 4: Vo{t)F A + Uo{t)G A = 0. (12) 



Due to the arbitrariness of u (t) and v (t), Eq. (10) implies g(x, t) = h(x, t). Ac- 
cording to the definition of h(x, t), one concludes that g(x, t) must be a real func- 
tion. 

According to Eq. (6), one can obtain uniquely ui(t) , vi(t) , U2(t) , V2(t) . Inserting 
them into Eq. (11), after some manipulations one has 

g (t)j t g 1 (t)-2g 1 (t)j t g (t)=0. (13) 
Solving for gi(t), we obtain 

9l (t) =C 1 (g (t)) 2 , (14) 
where C\ is an arbitrary constant. 

Inserting Eq. (14) into Eq. (11), one can obtain u 3 (t) when v 3 (t) is arbitrary. Re- 
placing ^i(£), Vi(t), u 2 {t), v 2 (t), u 3 (t) into Eq. (12), after some simplifications and 
let the coefficients of ip(t), ^-^{t), ^^{t) be zero (due to the arbitrariness of ip{t)), 
one has 

d = 0. (15) 

Thus Eq. (14) means that 

9i(t) = 0. (16) 

According to the definition of gi(t), gi(t) = -^g (x, t) \ x= -^(t), when ip(t) is arbi- 
trary, g(x, t) must be independent of x, so one has 

g(x,t)=g(t) = g (t), 
g 1 (t)=g 2 (t) = --- = 0. 

Finally, the conditions (15)-(17), together with compatibility condition (12) give 
the following equation 

4V (g (t)) 2 +(^g (t))g (t)-2 (j t g (t)) =0. (18) 
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From this equation we have 

e ±2Wot 

where A, B are arbitrary constants. For a proper choice of A and B, g (t) can 
change its sign at certain time (singularity point), which should be related to the 
Feshbach resonance. 

Thus, one can conclude that only when 

e ±2Wot 

g(x,t)=g(t) = Ae±AVVEt _ B , 

Eq. (2) can pass the WTC test. We will show below that this condition is sufficient 
for the integrability of Eq. (2). For convenience, it is also useful to introduce A = 
±2yfVo, i.e., Vq = ^A 2 . Thus g(x, t) becomes 

9 ^ t) = A^B- (19) 



It is very interest to note that under the present NLS equation, Eq. (19) is com- 
pletely satisfied the exact integrability condition proposed by Serkin, Hasegawa, 
and Belyaeva [ see Eq. (2) in Ref. (19)]. However, from the present analysis, Eq. 
(19) is a necessary condition for the Eq. (2) to pass the WTC test. A sufficient in- 
tegrability condition should be shown by further exploring a transformation which 
converts the GP equation (2) into a standard integrable NLS equation. 



3 Integrability 

Under the condition Eq. (19), Eq. (2) can be rewritten as 

d <9 2 2 g^e^ 1 2 1 

i—u (x, t )+Q^ u ( x > f ) + Ae 2\t _ B \ u ( x > *) I 2 u ( x > t) + ^ 2 x 2 u (z, *) = 0. (20) 



It should be noted that when A = 0, B = — 1 and A > 0, Eq. (20) is reduced to the 
case discussed in Ref. (11). 

Now we further prove that Eq. (20) is completely integrable. Our idea is to find a 
transformation to convert exactly Eq. (20) into the standard NLS equation 

i-^Q(X, T) + e^Q(X, T) + 5 \Q(X, T)\ 2 Q(X, T) = 0. (21) 
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In the literature, one often takes e = ±|, 8 = ±1 or e = ±1, 8 = ±2. 

In order to reduce Eq. (20) to the standard NLS equation (21), we look for a trans- 
formation of the form 

i>{x,t) = Q(X(x,t),T(t))e ia ^ t)+c{t) . (22) 

One also notes that other similarity transformations (23; 24) have been explored to 
reduce the nonautomonous NLS equation to the standard NLS equation. In com- 
parison to Eq. (22), the form of the similarity transformation in Ref. (23) is quite 
different from ours. The transformation parameters in Ref. (23) is independent of 
the space variables and the space dependence of the transformation is completely 
specialized by the explicit form of the similarity transformation introduced. In Ref. 
(24), another similarity transformation is introduced. However, this transformation 
reduces the nonautomonous NLS equation to a stationary NLS equation, as shown 
inEq. (2) in Ref . (24). 

After some tedious and technical calculations, one can obtain the explicit forms of 
the transformation 22. In the following we present them in different cases. 

Case I:A = 0,B = -1 

In this case, Eq. (20) becomes 

d d 2 1 

i— u (x, t) + — -u (x, t) + 2 g e xt \u (x, t)\ u (x, t) + -X 2 x 2 u (x, t) = 0, (23) 
at ox 2 4 

and the transformation reads 

T(t) = §f (e^-l), 

a(x, t) = —jlx £g ^ x H ^2 — J , 

Using this transformation, it is straightforward to check that Eq. (23) has the form 
of Eq. (21). Thus, according to a known solution Q (X,T) of Eq. (21), one can 
obtain a solution of Eq. (23) as follows 

<x,t) = V /i ^ eiAt Q( 2 f £eA *(-^ At )^( e2At - 1 )) e ^ 62X2 

(24) 

As noted in Ref. (25), this solution has some interesting features since it can take 
into account the fallout of particles form the BEC for A < and the amplification 
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of a soliton if A > 0. To show this, one can define a new variable Q'(x,t) = 
^2g e2 M u(x, t), which obeys the equation 



i^Q'x,t+-^Q'(x,t) + \Q , (x,t)\ 2 Q , (x,t) + ^X 2 x 2 Q , (x,t) = ^XQ'(x,t), (25) 

which includes an additional complex term at the right-hand side of Eq. (25). 

One notes that only when eSgo > 0, the above transformation is well-defined. In 
particular, for g > (g < 0), e = 1, 5 = 2 (S = -2) or e = ~, 8 = 1 (6 = -1) 
Eq. (23) denotes the standard NLS equation, which has been extensively studied in 
the literature. 

Remark: From the solution (24), it is obvious that if A > 0, the solutions of Eq. 
(23) can be compressed into very high local matter densities by increasing the ab- 
solute value of atomic scattering length. For the bright soliton (g > 0, e5 > 0) 
case, it has been discussed in details in Ref. (11). For the dark soliton (g < 0, 
s5 < 0) case, the situation is similar. In addition, if A < 0, the solutions can be 
expanded into very low local matter densities by decreasing the absolute value of 
atomic scattering length, as clearly seen from Eq. (25). 

Case II: A ^ 

In this case, the transformation can be obtained as follows 

T(t) - -2fiof£ ( i i_) 

1 W — XS 2 A \e 2Xt A-B A-Bj ' 
' . . \(Ae 2 ^+B)x 2 +^e^x+^- 

1 , / l , S(e 2Xt A-B)\ 

Likewise, this transformation can convert Eq. (20) into the form of Eq. (21). Note 
that only when 

5(e 2Xt A-B) 

— '- > 0, 

the above transformation is well-defined. The sign of e5 is determined by the signs 
of A and g , and the values of A and B. 



4 Explicit solutions of the GP equation 

As pointed out above, the solutions of the GP equation (20) can be obtained from 
the solutions of the standard NLS equation (21). This is a very efficient method 
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to find new solutions of the GP equation (20) since a number of solutions of the 
standard NLS equation (21) have been obtained by some classical methods such 
as the inverse scattering transform, the Lax pairs, the bilinear technique, and the 
Backlund transformation, and so on. 

In the following we first list some well-known solutions of the standard NLS equa- 
tion (21), which can be transformed into the solutions of the GP equation (20). As 
an example, we explicitly present the fundamental bright one-soliton solution of 
the GP equation. 

Fundamental bright one-soliton solution 

Q(X,T) = sech(X)ei T 

when g > 0, e — ^,5 — 1. According to Eq. (24), the fundamental bright one- 
soliton solution of the GP equation (22) can be written as 



u{x,t) = y/g^e* xt Q 

= Jg^e^ xt sech 



2 °" 2 



1* ^ x .2 4 9Q c \t x | f^O c 



For simplicity, below we only list the soliton solutions of the standard NLS equa- 
tion, the corresponding solutions of the GP equation (20) can be obtained by the 
same procedure shown above. 

Fundamental dark one-soliton solution 

Q{X,T) = q tanh(q X)e iq o T 
when q is an arbitrary constant and g < 0, e — — \,8 =1. 
Envelope solutions of the bright soliton 

Q(X,T) = hsechihX- hk 2 T) e -<- fc2X+ H fc 2 2 - fc i 2 ) T ). 
where k±, k 2 are arbitrary constants, and e = |, S = 1. 
Envelope solutions of the dark soliton 

Q(X,T) = k 1 tanh(k 1 X + k x k 2 T) e -<-^-|(fc2 2 +2fc 1 2 )T) _ 
where k±, k 2 are arbitrary constants, and e — — ~, 5 = 1. 

For bright and dark N-soliton solutions of Eq. (21), one can refer to Refs.(26; 27; 
28; 29; 30; 31; 32). 
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Besides the well-known soliton solutions listed above, in the following we discuss 
a new type exact solution of the standard NLS equation (21) by using the direct 
truncation method proposed in Ref. (33). These new solutions, of course, also give 
new solutions of the GP equation (20) 

Assume that 

Q(X, T) = rjihX + uxTytt 2 *^ 

be a solution of Eq.(21), where r)(.) is a real function. Some straightforward calcu- 
lations give 

u)\ = — 2 e k\k 2 , 

and 

- er, (0 h 2 - r) (0 uj 2 + e (J^r) h 2 + 5 ( V (£)) 3 = 0, (26) 
where £ = k\X + uj{T . By the similar procedure proposed in Ref. (33) we set 

(n _ a + bf(0 + cg(Q 

where a, b, c, p, q, r are constants to be determined together with ki, k 2 and lu 2 . The 
functions /(£) and g(£) are assumed to obey 

/(0 2 + s(0 2 = i. 

if(0 = g(0^-k 2 f(0 2 i ( 28 ) 



[|^) = -/(0 v /1 - fc W, 

Inserting (27) and (28) into Eq. (26), we get a polynomial of /(£) and Collect- 
ing the coefficients and let them be zero, one obtains a set of algebraic equations of 
a, b, c, p, q, r, k±, k 2 and uj 2 . Solving for the set of equations one can obtain many 
exact solutions. Here we only present three representative solutions. 

The solution I: 

Q(X,T) = ^-kk 1 cn(k 1 X-2ek 1 k 2 T, k) e -(-^ +£ ( fcl W-2^)r) ^ 

where ki,k 2 are arbitrary constants and e5 > 0. When k — > 1, we obtain the 
envelope solutions of bright soliton 

Q(X, T) = ^khsech (hX - 2 e k x k 2 T) e -i{~k 2 x+^-k^) T ) . 

When k — > 0, we obtain the period solutions 

Q(X, T) = \ (hX - 2 e k x k 2 T) e ~i-k 2 x+e(k^+k^)T) 

V o 



10 



The solution II: 



Q(X,T) = J-jkk 1 sn(k 1 X-2ek 1 k 2 T, k) e -<-^+^i W+^t)^ 

where ki,k 2 are arbitrary constants and eS < 0. When k — > 1, we obtain the 
envelope solutions of dark soliton 



Q(X, T) = y -jkhtanh (hX - 2 e k x k 2 T) e -i{-^ x +<^+^ 2 ) T ) 
When k — > 0, we obtain the period solutions 



Q(X,T) = J-jkk 1 sm(k 1 X-2ek 1 k 2 T) e -<- h ^+<^+^) T ) . 



The solution III: 



n(v r r\ I r fc snjkiX — 1ek x k 2 T, k) if 2k2X+£ ( 2k2 2_ kl * +2kl 2 k 2) T ) 
T) ~ V~25 kl 1 + cn (hX - 2 e k,k 2 T, k) 6 

where fci, k 2 are arbitrary constants. When k — > 1, we obtain the envelope solutions 

t\ r~^h tanh ( k ± x ~^ek x k 2 T) i(-2 fca x+ g (2fc a a +fci 2 )r) 



To our knowledge, the solution III is of a new type. 

Remark 1 Set k 2 = in the solutions I, II, and III, we get the stationary solutions. 

Remark 2 Obviously, when e5 > 0, the standard NLS equation has bright soliton 
solutions, and when e5 < 0, it has dark soliton solutions. 

Remark 3 As mentioned above, only when eSgo > 0, the transformation is well- 
defined. According to Remark 2 we know that, when g > 0, the GP equation (23) 
has bright soliton solutions and when g < 0, the GP equation (23) has dark soliton 
solutions. For the GP equation (20), one can use the Feshbach resonance to control 
the bright and the dark soliton in the BECs. This is because that if A 2 J°f_ B > 0, 

Eq. (20) is corresponding to Eq. ( 21) with e5 > 0, and if < °» Ec l- ( 20 ) 

is corresponding to the case of e5 < 0. The exact solutions of the GP equation 
with Feshbach resonance management presented in this paper have a significant 
contribution to the experimental study of the dynamics of the BECs. 

In conclusion, we discuss the integrability of the GP equation with Feshbach reso- 
nance management. By WTC test we find a condition under which the GP equation 
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is completely integrable. Meanwhile, we also find a transformation, which can con- 
vert the GP equation with Feshbach resonance management into the standard NLS 
equation. By this transformation, the well-known exact solutions of the standard 
NLS equation can be converted into the exact solutions of the GP equation. To our 
knowledge, this is the first time to give a correspondence from the exact solutions 
of the standard NLS equation to the exact solutions of the GP equation, and thus 
we can give abundant exact solutions for the GP equation. These solutions would 
be important for understanding the dynamical behavior of the BECs by using the 
Feshbach resonance technique. 
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